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Abstract 

We present results on the QCD equation of state, obtained with two different improved dy- 
namical staggered fermion actions and almost physical quark masses. Lattice cut-off effects are 
discussed in detail as results for three different lattice spacings are available now, i.e. results have 
been obtained on lattices with temporal extent of N-r = 4, 6 and 8. Furthermore we discuss the 
Taylor expansion approach to non-zero baryon chemical potential by means of an expansion of 
the pressure. We use the expansion coefficients to calculate various fluctuations and correlations 
among hadronic charges. We flnd that the correlations reproduce the qualitative behavior of the 
resonance gas model below Tc and start to agree with the free gas predictions for r>f. STc. 
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1. Introduction 

A detailed and comprehensive understanding of the thermodynamics of quarks and 
gluons, e.g. of the equation of state is most desirable and of particular importance for the 
phenomenology of relativistic heavy ion collisions. Lattice regularized QCD simulations at 
non-zero temperatures have been shown to be a very successful tool in analyzing the non- 
perturbative features of the quark-gluon plasma. Driven by both, the exponential growth 
of the computational power of recent super-computers as well as by drastic algorithmic 
improvements one is now able to simulate dynamical quarks and gluons on fine lattices 
with almost physical masses. 

2. Lattice parameter and scale setting 

We perform calculations on lattices of extent 16'^ x 4, 24^ x 6, using the p4fat3 ac- 
tion [P and 32'^ x 8 with both p4fat3 and asqtad actions. The latter calculations are 
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still preliminary and are currently performed by the HotQCD collaboration [2J. For the 
generation of gauge configurations we use the exact RHMC algorithm ^ . For each finite 
temperature calculation we perform a corresponding zero temperature calculation on a 
lattice of at least the size N'^, where N„ is the spatial extent of the finite temperature 
lattices. 

The simulations are done on a line of constant physics (LCP), i.e. the quark masses are 
kept constant in physical units. In practice, this has been obtained by tuning the bare 
quark masses such that the meson masses of e.g. pion, kaon and pseudo-scalar strange 
meson ss stay constant in the QCD vacuum as we change the value of the coupling. The 
strange quark mass was always fixed to its physical value, by fixing kaon and ss to their 
corresponding physical values [T]. We find that the LCP can, to a good approximation, 
be parameterized by a constant ratio of the bare quark masses. Most calculations are 
done on a LCP which corresponds to a pion mass of about 220 MeV. We do, however, 
also show preliminary results with a light quark mass corresponding to about 150 MeV 
for the lightest pseudo-scalar mass 4 . 

To set set the temperature scale in physical units, we determine two distance scales, 
rp and ri , from the zero temperature static quark potential 

\ / r=ro \ ^' / r=ri 

The ratio of both scales is only slightly quark mass dependent. It has been determined in 
both discretization schemes consistently, rp/ri — 1.4636(60) (p4fat3 Ij) and 1.474(7)(18) 
(asqtad [5]). The distance scales rg and ri have been related to properties of the cha- 
monium spectrum which allows to determine them in physical units. We use here ro = 
0.469(7) fm as determined in Ref. [6]. More details on the scale setting procedure, as well 
as the parameterization of the LCP are given in Ref. [l] . 



3. The equation of state 



Along the line of constant physics, at sufhciently large volume and at zero chemical 
potential, the temperature is the only intensive parameter that controls the thermody- 
namics. Consequently there exists only one independent bulk thermodynamic observable 
that needs to be calculated. All other quantities are than obtained by using standard 
thermodynamic relations. On the lattice, it is convenient to first calculate the trace 
anomaly in units of the fourth power of the temperature, Q^^/T*. It is easily obtained 
as a derivative of the pressure p/T*, with respect to the temperature. 

As the pressure is directly given by the partition function, p/T = \nZ, the calcu- 
lation of the trace anomaly requires only the evaluation of rather simple expectation 
values. According to Eq. [2] one then obtains the pressure by 
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Fig. 1. The trace anomaly, (e — 3p)/T* on Nt = 6 and 8 lattices. On the left panel results are obtained 
with the p4fat3 and asqtad actions, where Nr = 6 results are from [l] and [7] respectively, A'^t = 8 results 
are preliminary hotQCD results [2]. The right panel shows the low temperature part of the trace anomaly 
in detail, here we plot only results obtained by the p4fat3 action. Also shown on the right panel are 
quadratic fits to the data, as well as the trace anomaly obtained in the framework of the Resonance gas 
model. Light quark masses have been constrained to be one tenth of the strange quark mass (0.1ms), on 
the right panel we show, however, also preliminary results from simulations with physical quark masses 
(O.OSms) [1]. 

Here To is an arbitrary temperature value which usually is chosen in the low temperature 
regime where the pressure and other thermodynamic quantities are suppressed exponen- 
tially by Boltzmann factors corresponding to the lightest hadronic states; e.g. the pions. 
The energy density is then obtained by combining results for p/T'^ and (e — 3p)/T'^, 
respectively. 

In Fig. [T] (left) we show results for Q^^/T* obtained with the asqtad and p4fat3 
actions, respectively. The new N^- = 8 results [5] are compared to Nr = 6 results taken 
from [IJTJ- We note that the asqtad and p4fat3 formulations give results which are in good 
agreement with each other. In fact, in quite a large temperature regime the agreement 
for given lattice extent iV,- seems to be better than one could expect in view of the 
overall cut-off dependence that is visible when comparing results for Nt — 6 and Nr = 8 
more closely. They lead to a reduction of the peak height in O'^^/T'*, which is located at 
T ~ 200 MeV and to a shift of the rapidly rising part of in the transition region 

to smaller values of the temperature. 

In Fig. [T] (right) we show quadratic fits to the data obtained by the p4fat3 action, to 
highlight the cut-off effects. It is evident, that the Nr = 8 data is shifted relative to the 
Nr = 6 data by about 8 MeV at low temperatures, T ~ 160 MeV. This shift decreases 
to about 5 MeV at temperatures T ~ 190 MeV. 

Light quark masses have been constrained to be one tenth of the strange quark mass 
(O.lm^), on the right panel we show, however, also preliminary results from simulations 
with light quark masses of O.OStos. This again leads to a further shift in Q^^ of approx- 
imately 5 MeV at T ~ 190 MeV towards lower temperatures. 

We also compare the results for (e — 3p)/T'^ to results obtained from the hadron 
resonance gas model. Details on the resonance gas curve in Fig. [1] (right) will be given 
in [2] . The slope of (e — 3p) /T* obtained by the resonance gas model seems to be much 
smaller than the slope obtained by the quadratic fits to the data. Whether this points 
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Fig. 2. The ratio of pressure and energy density as well as the velocity of sound obtained in calculations 
with the p4fat3 action on A'^r = 6 [l] and Nt = 8 [2] lattices (left). The pressure in units of T* as funtion 
of temperature, for fixed and non-zero values of fig/T, obtained on Nt = 6 lattices (right). 



at deviations of the equation of state at lower temperatures from resonance gas behavior 
or is due to larger cut-off effects in the low temperature regime requires further studies. 
We note that the lattice spacing becomes larger at lower temperatures and violations 
of flavor symmetry, which are inherent to the staggered fermion formulations at finite 
lattice spacing, thus will become more important. 

The cut-off dependence observed in O'^^/T'* carries over to the calculation of pressure 
and energy density; the former is obtained by integrating over /T^ and the energy 
density is then obtained by combining results for p/T'^ and (e — 3p)/T^. This is apparent 
in Fig. [2] (left) where we show the ratio p/e obtained with the p4fat3 action on N^- = 6 
m and Nt = 8 [2^ lattices. Cut-off effects are still visible in the vicinity of the 'softest 
point' of the EoS, which is related to the peak position of (e — 3p)/T'^. We find that in 
the entire range of energy densities relevant for the expansion of dense matter created 
at RHIC, e < 10 GeV/inx^, the ratio p/e deviates significantly from the conformal, ideal 
gas value p/e = 1/3. 

This also is reflected in the behavior of the velocity of sound, — dp/de, which is 
shown in Fig. [2] (left) by dashed lines. It starts deviating significantly from the ideal gas 
value below e ~ 10 GeV/fm^ and reaches a value of about 0.1 in the transition region at 
energy densities e ~ 1 GeV/fm'^. Below the transition it slightly rises again, but note, 
that for very small temperatures c^, as well as p/e are sensitive the integration constant 
PoiTo) (see Eq.^. At present we have set po(To) = 0, for Tq = 100 MeV. 



4. Non-zero chemical potential 

At non-zero chemical potential, lattice QCD is harmed by the "sign-problem" , which 
makes direct lattice calculations with standard Monte Carlo techniques at non-zero den- 
sity practically impossible. However, for small values of the chemical potential, some 
methods have been successfully used to extract information on the dependence of ther- 
modynamic quantitites on the chemical potential. For an overview see, e.g. [Sj. 

We closely follow here the approach and notation used in Ref. J9]. We start with a 
Taylor expansion for the pressure in terms of the quark chemical potentials /iu,d,s, we 
obtain 
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The expansion coefficients c^f'i^{T) are computed on tlic lattice at zero cliemical poten- 
tial, using stochastic estimators. Some details are given in [lOlllj . We currently calculate 
the coefficients up to the 8th and 4th order, on A'^,- = 4 and 6 lattices, respectively. We 
find that cut-off effects are small and of similar magnitude as those found for the trace 
anomaly Q'^'^ [TT] . This was already anticipated by the analysis of the cut-off corrections 
in the free gas limit. Similar results for the asqtad action have been obtained in [T^. 

Alternatively to the quark chemical potentials one can introduce chemical potentials for 
the conserved quantities baryon number B, electric charge Q and strangeness S (/iB,Q,s), 
which are related to ^u,d,s via 

12 11 11 , . 

By means of these relations the coefficients cff^fj^ of the pressure expansion in terms of 
t^B,Q,s sltb easily obtained, in analogy to Eq. |4] 

In Fig. [2] (right) we show the pressure expansion up to the fourth order in /i^ /T. The 
results are obtained on Nr = 6 lattices. We find that corrections to the pressure arising 
from a non-zero chemical potential are dominated by the second order expansion coef- 
ficient for moderate chemical potentials of iib/T<2 and are of the order of 10-20% for 

T > Tr. 



5. Hadronic Fluctuations 

Quark number fluctuations are obtained from derivatives of the QCD partition func- 
tion with respect to the quark chemical potentials by the fluctuation-dissipation theorem. 
The Taylor expansion coefficients c^'ft , as defined in Eq. ffl can thus be directly inter- 
preted as quark number fluctuations at ^ = 0. However, quark fluctuations can not be 
detected directly in experiments due to confinement. Therefore we will consider fluctu- 
ations in terms of hadronic quantum numbers, i.e. baryon number B, electric charge Q 
and strangeness S, which are more easily obtained by experiment. These fluctuations are 
related to the Taylor expansion coefficients c^fl^^, as given in Eq. [6l A recent overview 
on the physics of fluctuations in the context of heavy ion collisions was given in Ref. [13j . 
In general, the quadratic fluctuations X2 zero chemical potentials can be obtained 
from the second order coefficient c? 



X 



^2 

1 S^lnZ 



2VT^ ^{^lx|Tf 



Mb,Q,S— U 



mxY), , (7) 



where 5N = N — (N) denotes the normalized net-density, and (. . .)o indicates that the 
expectation value has been taken at fiB,Q,s = 0. Under such conditions, baryon number. 
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Fig. 3. Quadratic and quartic fluctuations of baryon number, electric charge and strangeness, normalized 
by their corresponding Stefan-Boltzmann value. The results on N-p = 4 lattices (open symbols) and 
Nr = 6 lattices (full symbols) are in good agreements. 

electric charge and strangeness vanish, and we have SN ~ N. We define quadratic and 
quartic charge fluctuations by 



X2 



1 



1 



((iVi)o-3(7Vi)^)^24cf , 



respectively, and correlations among two conserved charges by 
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where X,Y e {B,Q,S}. 

In Fig. [3] we show results for quadratic and quartic fluctuations of B, Q and S. The 
quadratic fluctuations X2 '^'^ I'ise rapidly in the transition region where the quartic 
fluctuations Xa''^'^ show a peak. The peak height is more pronounced for the baryon 
number fluctuations than for fluctuations of the strange quarks. 

We compare the results obtained on lattices with temporal extent N^- = 4 and 6. 
We notice that they are in general compatible with each other, especially in the high 
temperature phase, where both quadratic and quartic fluctuations approach the Stephan- 
Boltzmann limit quickly. The transition temperature has been previously determined 
to be Tc = 202 MeV and 196 MeV on = 4 and 6 lattices respectively [H]. We 
thus conclude that at temperatures of about 1.5Tc and higher, quadratic and quartic 
fluctuations of i?, Q and S are well described by the ideal massless quark gas. 

At low temperature, hadrons are the relevant degrees of freedom. The hadron resonance 
gas (HRG) model has been shown to provide a good description of thermal conditions at 
freeze-out. We thus compare the fluctuations in the low temperature phase with a HRG 
model, where we include all mesons and baryons with masses smaller than 2.5 GeV from 
the particle data book. 

In Fig. m we show the ratio of quartic and quadratic fluctuations for B, S and Q. In 
the HRG model, xf /xf easily obtained in the Boltzmann approximation, which is 
valid for a dilute baryonic gas in the temperature range of interest. One finds that all 
details on the hadron mass spectrum and temperature dependence cancel and the result 
is a constant, given by the unit square of the baryonic charge (one for all baryons). This 
is in fact reproduced by the lattice results shown in Fig. |4] (top left). 
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Fig. 4. Ratio of quartic and quadratic fluctuations of baryon number (B), strangeness (5) and electric 
charge (Q). The two curves of the HRG model in the top right figure correspond to charge fluctuations 
with physical pions (upper) and inflnitely heavy pions (lower curve). 
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Fig. 5. Pairwise correlatios of the conserved charges baryon number (B), electric charge (Q) and 
strangeness (S) as function of the temperature, normalized to the quadratic fluctuations of B and Q 
respectively. 



The ratio of quartic and quadratic fluctuations for S and Q are more complicated even 
in the Boltzmann Umit, since hadrons with different electric/strange charge give rise to 
different contributions to the corresponding fluctuations. For strangeness fluctuations, 
shown in Fig. |4] (bottom), the Boltzmann limit is still a good approximation; but for 
electric charge fluctuations, the pion mass plays an important role. In order to check for 
the sensitivity of electric charge fluctuations on the pion mass, we show in Fig. |4] (top 
right) results of a HRG model calculation with physical pion masses and without the 
pion sector, i.e. for infinitely heavy pions. 

In Fig. [HI we show the various correlations Xi? j Xn^ ^-iid Xii normalized to quadratic 
fluctuations X2 ^'^^ X2 respectively. The results from Nr = 4 and 6 lattices agree 
with each other very well, and they are compared with the HRG model in the low 
temperature phase and Stefan-Boltzmann limit in the high temperature phase. We flnd 
that the correlations reproduce the qualitative behaviour of the HRG model below 
and again start to agree with the free gas predictions for r>1.5rc. The rapid suppression 
of fluctuations above Tc as well as the agreement of the correlations with the free gas 
predictions in the high temperature phase, suggests that in the quark-gluon-plasma the 
baryon number and electric charge are predominantly carried by quasi-particles, with the 
quantum numbers of quarks T5]. 
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6. Conclusions 

We have presented results on the equation of state on lattices of Nr = 4, 6 |1|7| and 8 
[2] obtained with two different kinds of improved staggered fermions. Hadron masses have 
been kept constant in physical units and are chosen such that we have a physical strange 
quark mass (jUs) and 2 light quarks with a mass of mi = O.lmg. We also presented some 
preliminary results with smaller quark masses mi = 0.05ms. We find that our two actions 
lead to a consistent picture of the thermodynamics of QCD and find in particular for 
the Nr = 8 results only small cut-off effects. We have calculated the equation of state 
as well as the velocity of sound and find the softest point of the equation of state to be 
{p/e)min — 0.09 at energy densities of 1 — 2 GeV/fm'^. 

Furthermore, we calculated corrections to the equation of state arising from a non-zero 
baryon chemical potential, by means of a Taylor expansion of the pressure. Using the 
expansion coefficients, we have analyzed the quadratic and quartic fluctuations of baryon 
number, electric charge and strangeness, as well their ratios. We find these quantities to 
be in good agreement with the free gas results at temperatures of T > 1.5Tc. Below 
Tc, qualitative features of the resonance gas are reproduced. The ratio for the baryon 
number is closely related to the second approximation of the convergence radius of the 
Taylor series of the pressure with respect to the baryon chemical potential. 
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